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The variation of the perturbative 3-curvature parameter, £, is investigated 
. in the period of reheating after inflation. The two-field model used has the 

inflaton, with an extra scalar field coupled to it, and non-linear effects of both 
fields are included as well as a slow decay mechanism into the hydrodynamic 
04 ' fluid of the radiation era. Changes in £ occur and persist into the succeeding 

cosmic eras to influence the generation of the cosmic microwave background 
fluctuations. With appropriate coupling the extra field rises from almost zero 
at the end of inflation to parametric resonance which is important in giving 
i-£h ' ultimate changes to £. Inflaton particles are always produced in important 

quantities and their interaction with the extra field may limit its paramet- 
ric resonance. When the field coupling is such that parametric resonance is 
,^Zh • small or absent then non-linear effects of the produced inflaton particles by 

themselves also give ultimate changes in £, though these are small. 

• i— i - 

X 

I. INTRODUCTION 

In the cosmic inflationary scenario parametric excitation of extra fields may mediate the 
transition of the inflaton field into a cosmic radiation era of a gas of relativistic particles |T]-|| . 
There has been discussion of the possible role of non-linear effects in such transitions flOHl4 



depending on their magnitude such effects may be important in connecting the quantum 
excitations of the inflationary era inflaton field to the cosmic microwave background radiation 
fluctuations (CMBRF). When non-linear effects are not important this connection has often 
been made rather directly, for example by using the time-invariant properties of a parameter 



C, associated with spatial curvature perturbations ||15|-[L7j1, for wavelengths appropriate to 
the CMBRF. To study such non-linear reheating effects on the CMBRF we calculate using 
the simplest theory as a test-bed with just one additional scalar field subject to parametric 
resonance [pT0| ,|lT[] . 

Thus we have a model with an inflaton field ip = ipo + ipi, where ipo denotes the 'classical' 
part of the inflaton field and tp% the part arising from primordial quantum perturbations 
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of if. We use the perturbed FRW metric in longitudinal gauge which appears with the 
perturbing fields if and $ as 



ds 2 = a(r) 2 (l + 2$)dr 2 - a(r) 2 (l - 2ip)S ij dx i dx j . (1) 

In the usual case of absence of space-space off-diagonal elements of the stress-energy tensor 
if and $ are equal. The one additional field x is taken to have the well-used interaction 
term \g 2 f 2 x 2 EHH- Persuasive arguments have been given []12-14]| that x is strongly sup- 



pressed before reheating begins. It is then reasonable to treat x as a field arising from its 
own quantum fluctuations without any primordial classical part. We thus have 3 initial 
quantum fields ip, if\ and %. The equations of motion imply that in the canonical decom- 
position of quantum fields into momentum component creation and annihilation operators, 
the operators for if and if\ are the same and those for x are independent. 

As shown by Polarski and Starobinsky [Tj|] the quantum fields decohere into classical 
stochastic fields which nevertheless keep important properties of the initial quantum fields; 
for instance the ensemble averaging of the products of stochastic fields give similar results 
to the vacuum expectation value of the corresponding quantum fields. Since we cannot real- 
istically invoke a non-perturbative and exact quantum mechanical treatment such averages 
necessarily arise in our treatment of the equations of motion; for example from \g 2 f 2 X 2 a 
term of the form g 2 ifoX 2 appears in the equation of motion of the classical homogeneous 
part <fo of the inflaton requiring for consistency the introduction of the ensemble average, 
< X 2 >, of X 2 - 

This feature is important, for example in the case of the homogeneous Einstein equations 
where the contributions due to if\ and x during the reheating period may become the 
dominant ones. Also in the non-homogneous Einstein equations if \ cannot be treated as 
a perturbation but should be viewed as the non-homogeneous part of the field if. On the 
other hand the field ip(x, t) will be treated as a perturbation and in our numerical simulations 
checks will always be made to ensure that ip(x,t) is always much smaller than unity. The 
justification of this procedure has been discussed previously [flQ| , p~T|] . 

Our model also incorporates true heating through an elementary decay mechanism, with 
decay constants and T x , of the scalar fields into a hydrodynamical relativistic (radiation) 
fluid. The system thermalizes well with if becoming approximately constant as expected 
from analytic calculations. We calculate, using well-known procedures, the cosmic microwave 
background fluctuations starting from the values of if or ( immediately after thermalization 
thus incorporating any changes that may have occured in the transition from inflation to 
the radiation era. 

We shall see that we have to deal with two different regimes depending on the relative 
values of the particle or field energy densities attained in the transition. On the one hand we 
have the classical inflaton field, fo, whose density declines plus the inflaton particles field, 
if i, whose density increases; these densities become close to each other and stay so until they 
thermalize. On the other hand we have the x density which at some time in the transition 
increases rapidly from almost zero to some non-negligable value. One regime, A, is that 
where the x density stays below the other densities; the second regime, B, which is induced 
by larger couplings than those of A is where the the x density approaches rivalry with the 
others. Both regimes can be induced by a large number of different values of the parameters 
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g 2 , the coupling strength, and of and T x which vary the rapidity of thermalization. In 
both regimes ( changes but in A, where the x density remains comparatively modest, the 
change is constant and smaller while in B the change factor varies. 

We know from analytic argument that when there are non-linear effects of the fields the 
usual theorems showing ( to be constant fail, for example as shown in section p .f"E| below. 
Because of the complexity of the equations we do need numerical computation to verify and 
quantify the change in ( and any importance it may have for calculation of the CMBRF. 
This paper attempts to do so for a particular model of preheating and reheating. 

The paper has been organised as follows. In the next section we precisely define the 
model for the reheating period and give the resulting equations of motion. The initial 
conditions arising from the inflationary period are also given. We call attention to the fact 
that despite the initial extreme suppression and smallness of x this fi e ld still develops 

an important resonating behaviour. The theory can be formulated in different ways because 
of the redundancy of the many equations of motion and some small such changes have been 
made with respect to our previous papers [|K||ll[] with the purpose of increasing the accuracy 
of the numerical computations. In section III we discuss the results of the computations and 
in section IV we make the connection with the cosmic microwave background fluctuations. 



II. THE REHEATING EQUATIONS 
A. The Equations of Motion 

The scalar field part of the model is specified by the Lagrangian 

L = J d 4 x^g-{^> a v, a + \x a X^ - V(<p) - V( X ) ~ V int {cp, x)} (2) 
where ip = <po + tp\.. We take the field potentials during reheating to be 

v(tp) = ~mV; V( X ) = \m\ 2 - vu<p, x) = \gW- (3) 

The perturbed longitudinal gauge metric is given by Eq.(fj]), the dynamics of the evolution 
of the scale a(r) by the Einstein homogeneous equation (|T2|) below. It is characteristic of 
our treatment that we always include the metric perturbations. As in our physical system 
there are no space-space off-diagonal elements of the stress-energy tensor R, $ = if> and we 
shall denote the metric perturbation field by if). A friction mechanism will also be adopted 
to simulate the decay to the thermal gas of relativistic particles dominating subsequent 
radiation era. These terms take the simple form T^cf)' 2 and F x x' 2 and they can be consistently 



x The ensemble averaging of non linear terms, described below and more fully in ref. suppresses 
any such terms which might otherwise be present in our non-linear theory. Corrections to the result 
$ = ip would be present in a full quantum mechanical treatment, but as stated in the introduction 
that is beyond the scope of the present work. Finelli and Khlebnikov [19| also find such corrections 
in a different context where the fields are classical rather than in a stochastic ensemble as here. 
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introduced into the equations of motion in such a way that the Bianchi identities are satisfied 
0- 

The stress-energy tensor and the equations of motion involve the averages of bilinear 
products of stochastic fields. Using as an example <pi these are defined as follows: 

(Vi(x,T) 2 ) = (2n)- 3 Jd 3 k VkV t (4) 

where (p^ are the k-mode functions entering in the expression 

— — r[e(k)^ k (r) exp(zk.x) + e*(k)<p£(r) exp(-zk.x)], (5) 
(27r) 2 

with e(k) and e*(k) being time-independent 5-correlated Gaussian variables (TB| such that, 
where (....) denotes the average, 

<e(k)e*(k')> = ^ 3 (k-k'); (6) 

with averages of ee and e*e* products being zero. The field ^(x, r) is expressed as Eq.([5D 
with the same gaussian variables but the mode functions being ip^r). The stochastic gaus- 
sian variables of x are independent. Thus, for example, {<pi{x, t)x{x, t)) is zero while 
{(fii(x, r)ip(x, r)) and (x(x, r) 2 ) are non-zero flTOl . 

The thermal gas of relativistic particles is described by the hydrodynamical variables 
p, the density, and p, the pressure, related by the equation of state p = |p. As with the 
inflaton field the density has two components 

p(x,r) =po(r) +pi(x,r), (7) 

being the homogeneous background (classical) density and the non-homogeneous quantum 
or stochastic part. As will be seen below the equations of motion imply that the stochastic 
variables of p\ are the same as those of <pi and if). We assume that p is zero at the beggining 
of reheating and only develops subsequently through the friction mechanism T^' 2 and T x x /2 . 
In terms of consistency of the friction model the assumption that po only becomes significant 
during the reheating is not unreasonable because the friction only becomes significant when 
the fields are oscillating rapidly. 

From the Lagrangian and from the hydrodynamical variables p and p we define the 
energy-momentum tensor 

T* = + X*X,v ' [\v a V,« + \x a X, a ~ V(cp) - V( X ) - Vim - PK + (P + V)u»u v (8) 

where is the 4-velocity of the fluid. 

The homogeneous part of the total density, Pt(x, t), is defined as Pot(t) =< T® >, 
giving, after ensemble averaging, 

Pot(t) = 7^iv + a'm^l + (^)) + a 2 M 2 ( X 2 ) + { V { .} + ( X 2 > - Vo^)] + Po, (9) 
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(10) 



m 2 =m 2 + g 2 ( X 2 }. (11) 

The definition of por is slightly more subtle. Although we shall be using the Friedmann 
equation to calculate a(r): 

(a'/a) 2 = ^fa 2 p 0T (r). (12) 

we should ensure that this equation is indeed compatible with the second order (time-time) 
Einstein equation for a" /a which is: 

a" /a - (a'/a) 2 = -^a 2 (p or (r) + 3p QT (r)). (13) 



Deriving Eq. ([r2j) and comparing with Eq. ([T3|) we find 

1 



PoAr) = ^[v- - a 2 m% 2 + (tf)) - a 2 M 2 ( X 2 ) - V <<^>] + 

4 



2a 2 ' 

\-iro„-2,„/ 



p /3 - (a'/a)- 1 ^- 2 ^^') + ~{p^')\. (14) 

V- = ^ + <¥?> + <X' 2 > - <tf,«>/3 - (xJ)/3. (15) 
This differs, in the addition of the last term in square brackets, from the expression obtained 



by using Pot(t)5j = —(Tj) which we have used previously fT0|l , but it has the additional 
virtue of precisely maintaining the continuity (or 'energy conservation') equation. However 
in the equations of motion the expression for p T only occurs in the metric perturbation 
equation (Eq.(^) below) and in fact the addition of the last term seems to make little 
difference to the numerical results. 

Besides the Friedmann equation there are equations of motion for the other homogeneous 
variables^], these being y? and po : 

<p% + 2(a'/a)<p' + a 2 m 2 <p - 4(VVi) - 4(^VVi> + 2a 2 m 2 (^!> = ~<^M + 2(M» ( 16 ) 

p' + A(a'/a) Po = a-'T^' 2 + (<p*)) + a~ l Y x { X n ) - 2a~ 2 ^V 2 Vi) + M^'pi) (17) 

There remain the spatially non-homogeneous equations which we write in the A;-component 
form. These components are specified as the complex mode functions Xki ¥>k, 4>k and pu-. 
Their wave number dependence, given by the succeeding equations, is only on k = |k|. 
There are 4 equations which we take as those for Xk, <£k and the time-time and space-space 



'Both equations can be obtained from the Bianchi identities as shown in Appendix A of ref. [1C] 
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Einstein equations for ip^. (These being satisfied the time derivative equation equation for 
p lk = 3pik is implicit through the Bianchi identities.) 

Xl + 2(a'/a) X ' k + (k 2 + a 2 M 2 ) Xk = -aT xX ' k (18) 

where M 2 is a function of r given by 

M 2 (r) = M 2 + g 2 Mrf + <?%i(x, r) 2 ), (19) 

tp'l + 2(a'/a)<p' k + (k 2 + a 2 m 2 )(p k - ^'^' k + 2a 2 rh 2 ip ^ k = -aT v (ip k + 2ip' il> k ), (20) 



the equations actually holding for each separate value of k; thus consistency of Eq . (f20|) 
justifies our previous statements that ifo., the mode functions of the metric perturbation, 
should be associated with the the same stochastic variables (or quantum operators) as ip^. 
The space-space Einstein equation for the metric perturbation can be written 

+ 3if/ k (a'/a) + /a) 2 = AnGa 2 p T k + ^-a 2 (p 0T + p 0T )^k (21) 
where pxk is the /c-component of the non-homogeneous part of the total momentum: 

PTk = a~ 2 [-T] + ip k + (p' (p' k - 2(p' k (if)(p' 1 ) - a 2 fn 2 (pQ(p k ] + p lk , (22) 

^+ = ^ + (<P?) + (X' 2 ) + (^)/3 + (x3)A (23) 
The time-time Einstein equation is 

3ip' k (a'/a) + i) k {k 2 + 3(a'/a) 2 ) = 4irGa 2 p Tk . (24) 

PTk = a' 2 [-(r] - (<pi,i<pi,i) - (x,iX,i))^k + f'of'k ~ 2^' k {^^'i) + a 2 rh 2 (p (f k } + p lk . (25) 



Rather than using the mode equation for p' lk to find p\ k we find it from Eqs. (24.25) and 
then use the equation of state p\ k = \pik to solve the ip equation ( |2"1~| , |2"2"| ) . ^ 

In addition there is the time-space Einstein equation which is not independent of the 
other two Einstein equations but which acts as an equation of constraint on the initial values 
where we use it to fix the value of ip' at the beginning of reheating. It takes the form 

i)' k + (a'/a)^ k = 4vrGVo<^ (26) 



3 As are ip k and ifi, x is complex and the real and imaginary parts are not constant multiples of 
each other for smaller occupation numbers of a mode 

4 The pi k mode equation is still maintained; the Bianchi identities hold and imply the p\ k mode 
equation when the other equations of motion are explicitly fulfilled, as they are |1C| ]. 
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B. Initial conditions 



We need to have the values of H, the fields and their time-derivatives at the end of 
inflation to supply the initial conditions for the reheating equations of motion . We have 
chosen to use a specific inflationary model having the advantage that the solutions are 
analytically expressible. This is power-law inflation with an exponential potential: V = 
U exp(—\(p) where U, A are constants. This potential of the inflation era stands in place of 
the potential V(<p) = \w?Lp 2 of the reheating era but otherwise the Lagrangians are the same 
and in particular both have the interaction potential V int ((p, x) — \g 2{ P 2 X 2 - We have ensured 



the correct continuity of the potential by imposing the Lichnerowicz conditions |2t| , as well 
as using the equations of motion and constraint appropriate to each era at the boundary 



That the x-field be strongly suppressed in the inflation referred to previously 

12| ^14| , is an essential feature in maintaining the analytic form of power-law inflation. The 



many efold increase of a during the inflationary era indicates an exceedingly small value for 
Xk at the beginning of reheat 0^0' which may be as little as 10~ 50 of its initial value. 
This is an important qualitative feature of our initial conditions. We have taken the initial 

— 1 /2 

value of x to be 10~ n m P ^ ncA , where n is of the order of 30. (For such small mode functions 
in the beginning of reheat the transition to classical stochastic functions cannot yet be made 



and the complex formalism, which is applicable also to the quantum case []T8| , [10f| should be 
retained. 

For power-law inflation the scale factor a oc (t^ — r) p and we are free to choose p by 
specifying A in the potential flOfl . In the numerical results quoted we have chosen p = —1.1 
(p = —1 corresponds to exponential inflation). 

Knowing the analytic solutions of the unperturbed equations enables the deduction of 
ip k and (fik in the inflation era starting from the initial quantum oscillations ||10|| . These are 



ipk = a l J^^O(^\j\(p,' k -a^ k )/k 2 (27) 



<p k = a~ 1 (2k)- 1 / 2 [fi k + a 7 04 - ap k )/k 2 ] (28) 
where 7=1 — a'/ a 2 = (p + l)/p and 

Hi + (k 2 - a"/a)p k = (29) 
Since a" /a = ^-(t, — r)~ 2 , Eq.Q2Tf) can be expressed as a Bessels equation in terms of 

y = k{ n - r) = \p\ (30) 
an 

and the solutions which satisfy the asymptotic conditions for a scalar field are 

I Try 1 1 

= y -!}{ J v{y) - iY v(y)) exp[-i(-i/7r + -tt)] (31) 

where v — 1/2 — p. 
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For k/a <C H (many times over fulfilled by the k relevant to the CMBRF observations) 
/ifc(y) oc y p oc k p and if this fc-dependence were to pass unaltered through the reheat era into 
the radiation era it would yield a power-law spectrum scaling as k n ~ x where n, the spectral 
index, is given by n = 2p + 3. 

From the conditions of continuity it follows that at the junction of inflation and reheat 

H = V^rGrwp y/p/(2p-l) (32) 
and this, using Eq. (|30|) , determines the value of y at the junction 



III. SOLVING THE REHEATING EQUATIONS 

One of the main aims of the present work is to see whether, and if so how far, parametric 
resonance of a second scalar influences the fluctuations of the CMBR. 

An important parameter is ( k HHim wm ch, for k 2 /a 2 <C H 2 , is equivalent to the curva- 



ture perturbation lZ k - ( is defined in terms of the gauge- invariant metric perturbation $, 
equal in our case to the longitudinal gauge perturbation ip; in terms of this variable we have 



C k = ^(H- 1 ^ k + ip k )/(l + w)+^ k ^-TZ k (33) 

where w = Pot/ Pot, the total (homogeneous) pressure to density ratio. If k 2 /a 2 <C H 2 
then ( is constant during the radiation era and constant if there be only one scalar field, 
the inflaton, during inflation and also during reheating if the inflaton perturbation can be 
treated to first order. That is, as it is often expressed, it is constant from horizon exit 
during inflation to horizon re-entry around the beginning of the matter era. The question 
arises whether, with another scalar field having parametric resonance during reheating, can 
( k exit reheating with a different value from that with which it enters [ |l2| , f23"| , p4f As we 



shall see, for values of g/m bigger than about 5000 there is parametric resonance, with \ 
comparable to the other scalar densities, and ( is larger than its initial valueQ and also there 
can be contributions to this effect from higher order inflaton perturbations. Then, knowing 
( at the end of reheating (that is at the beginning of the radiation era, asessed as in the 
succeeding paragraph) we can use known methods to find some properties of the fluctuations 
and assess consequences of the change in £ in reheat. 



Another method pl| , p2| is to trace the metric perturbation ip k , for the appropriate very 
small values of k, from the end of reheating to the matter era where the Sachs- Wolf effect 
comes into play. Once the system thermalizes, with the gas of relativistic particles obeying 
an equation of state of the form p = p/3, we can check that as the reheating era tends to 
the beginning of the radiation era then ip k = bl + b2/a 3 ,a(r) « r, is a solution of Eqs. 
(PUMPI), f° r k 2 ~ 0; where the constant bl has to be determined from the computations. 



This is a correction to some numerical values given in our previous work. Due to a misprint in 
the computing code some decreasing values of ( for large couplings were, wrongly, found. 
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The attainment of constant value for ipk^o and pa A will be taken as indicating the beginning 
of the radiation era; ip can then be traced throughout the radiation era from its analytic 
form which holds therein. 

In this paper we evaluate the power spectrum of the fluctuations and find the quadrupole 
moment (which can give the COBE normalization) using the first method above. 



In contrast to our results given below, the work of refs. |L2| finds very small effects of £ 
variation by expressing ( in terms of the non-adiabatic pressure variation. This comes using 
the reasonable physical assumption - based on the approximate validity of the Robertson- 
Walker metric - that after smoothing on a cosmological scale below the ones of interest the 
spatially inhomogeneous variables such as pk,Pk can be treated as perturbations. There 
is no evident agreement between that expression for ( and Eq. (|34D below which uses no 
assumptions on pk,Pk being perturbative. However, as we discuss in sections I and V, in our 
work we do consider the variation ip from the Robert son- Walker metric to be perturbative. 

A. Small coupling constant, g/m 



Having derived in Sec. [II B| the expressions allowing the magnitude of the relevant modes 
at the beginning of reheat to be calculated, thus fixing the input of ip, tp\ and x t° the reheat 
era, we may proceed with the numerical integration of the reheating equations of Sec.|ll A|, 



taking the integration up to a point where we may consider the system to be well thermalized. 
The output of this numerical development will then be used in Sec.[TV|. 

Unless otherwise stated, the results are for p = —1.1, M/m = 0.02, </> = .3mpi anc k at 
the beginning of reheat and the inflaton mass m = 10~ mpi anc k [m P i anck = G- 1 ' 2 ). The 
values for the coupling constant g, in |g 2 y? 2 X 2 > wm De quoted in terms of g/m. 

The frictional decay constants T x , have been investigated in the range 0.1 > T x /m > 
and 0.005 > T v /m > 0, larger values tending to suppress the resonance. However the results 
presented here will be for small values of the order of 10 -4 . The reason is that then the 
parametric resonance has had time to reach its full effect before the (po, (pi and x decay. For 
larger decay constants the transition to the radiation era occurs when ( is still fluctuating 
very strongly leading to more random ( values in the radiation era. It is significant that 
without this physical assumption the theory would be very randomly dependent on the 
precise magnitude of the coupling. 

In this subsection we consider values of g/m < 4000 and some typical results are given 
in Table I. The main conclusion seems to be the fact that the end of reheating Q values 
of our variables ipk^o an d (k^o are practically independent of g/m (though not of the T's 
which govern the rate of thermalization) . This was checked varying g/m including g/m = 0, 
obviously a no- resonance case, up to nearly g/m = 4000. Such results, at first surprising, 
may be understood looking carefully at the energies associated with the fields ipo, y?iand x- 

The graphs in Figs, 1,2 and 3 show the reheating transition for g/m = 3000 and T x = 

= 5 x 10~ 5 . Fig.l shows the energy densities of the inflaton particles (pi and of Xi 
displaying the parametric resonance. The scale is such that in all the graphs not all the 



3 The suffices in and out denote the beginning and end of reheat respectively 
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rapid oscillations are adequately displayed. Fig. 2 shows the considerable ( oscillations 
during the transition and that it is rather constant in the initial period when, as we can see 
from Fig. 3, the density of ipi is small compared to <^ . Despite considerable oscillations to 
larger ( the final result for these small decay constants is that ( decreases and by relatively 
little. 

Although < x 2 > increases by many powers of ten during the resonance (more than 
10 20 , for g/m = 3000, m = 10~ 7 , the energy associated with this field is, nevertheless, 
always smaller than the energy of the fields (fo, <pi, which thus control the dynamics of the 
system. The most striking feature is that the case g = having essentially zero \ field 



produces the same change in This has already been found and discussed in ref. ||TT| for 
the case where there is no particle decay (r x = = 0). There the analytic expression for 
the rate of change ( is non-zero due to the binary (non-perturbative) terms in the inflaton 
particle field, tpx even if the x fi e ld were zero: 

-if (1 + w)( = # + Hip + 2H^ - $ + Hip)—\n(p h + p h ), (34) 



/ , x _ Wo + (Wi) + (xx) + + <X,»X,i)/3 

dt iU[Ph +Ph) - Z vl + (Vl) + (X 2 ) + (^ M >/3 + (x^,i>/3 ' [ } 

Thus the presence of the inflaton particles, comparable in their contribution to the energy 
to the residual homogeneous inflaton field, generates a change in ( in the reheating - or more 
exactly the preheating - period. However as we see from Table I, for the small value of 
chosen (for the reasons given above), the ultimate change in £ is quite small. 



B. Large coupling constant, g/m 



We see in Table II that ( increases in this regime: ( out > Q n . The graphs in Figs, 4, 
5 and 6 show the reheating transition for g/m = 4700 and T x = = 5 x 10~ 5 . We may 
note that: firstly the parametric resonance in x is greater as one would expect with a larger 
coupling; secondly that ( out now varies with g. Thirdly the fluctuations in ( tend to be 
larger, with a considerable increase in a smoothed value of ( before true heating takes over, 
the increase taking off when the x density attains the (p± density as can be seen from Figs. 4 
and 5. Again we should emphasize the role of the decay widths; if these had been much 
larger in that case then the final ( would have been much changed. 



C. The sudden approximation 



It may be of interest to compare our results with those of the sudden approximation, 
where there is no use of the reheating mechanism, especially because of the independence of 
g in a range of smaller couplings set out in section [111 A| . The matching conditions, deduced 
by Deruelle and Mukhanov |20| from the Lichnerowicz conditions, imply |21] that ip an d the 
function 
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r = (fa' 1 + if; - \v 2 ipa~ 2 )/{l - a or 2 ), (36) 

as well as a = a' /a, be continuous on the sudden transition hyper surf ace, say r = r 2 . At the 
end of inflation, that is on one side, denoted r 2 _, of the transition hypersurface, the values 
of if> an d r may be written 

™ = ^TI E ' (38) 



where S oc y p is given by the inflation model as in section IIB. On the other side of the 



hypersurface in the radiation era, we may write the general solution for if) 

ipk = [Bi(e cos e — sine) + B 2 (cose + e sine)] /ae 2 , (39) 

where e = e — e 2 , e = and e 2 = /cr 2 /v^3 is the value of e at the beginning of 

the radiation era. (We take the convention on r that a = cr, c constant, in the radiation 
era.) From Eq. fl3"9]) we can obtain the values of ip and T on the radiation era side of the 
hypersurface denoted r 2 +. Then solving 

^(r 2+ )=^(r 2 _), (40) 
T(r 2+ ) = r(r 2 _), (41) 



gives 



2 . P+l 



5ie 2 = -2a 2 S + 2e^a 2 ^— — S, (42) 

2p + 1 

fi 2 = 2a 2 S - e^a.^^S (43) 

For values of k relevant to the CMBRF we can expand if) in powers of e through the pure 
radiation era: 

i,(r)=f a + 0{e 2 2 /e 2 ) (44) 

V*. = = ?S (45) 

3 e 2 a 3 

where in the last equation we have used the fact that both a and e are proportional to 
r. Thus, in the radiation era, if) changes rapidly at the beginning and after a few efolds 
of a takes up a constant quasi- asymptotic value, ip qa . This constancy of if> in most of the 
radiation era is a general property and not peculiar to its initiation by a sudden transition. 
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For the sudden transition one finds that if) qa is considerably bigger than if) in , the value at 
the end of inflation: 

Ipqa = o P+ 1 ^in = 7-33?/W (46) 

3 p+1 

The corresponding quasi-asymptotic value found in our gradual reheating calculations for 
g/m < 4500 is surprisingly close: 

VV « 7.0^ n - (47) 

A similar explicit evaluation of ( in the sudden transition case gives from the beginning 
of and all through the radiation era 

C=--S = S = 0„- (48) 
e 2 a 

In the gradual reheating calculations for g/m < 4500 

C = taut = -88Cin- (49) 

Thus when parametric resonance is small our gradual transition calculation and the sudden 
transition result approximate each other to within about 12% for ( and 5% for if>. 



IV. CMBR FLUCTUATIONS; QUADRUPOLE MOMENT 

We consider that the ( effects we have calculated for a non-linear theory would apply, 
qualitatively, to a large range of parameters. Nevertheless we should show that the param- 
eters used here are quite near to ones which correctly normalise the theory to the cosmic 
microwave background fluctuations as required for modelling the real world. 

So in this section we shall be concerned with the calculation of the scalar perturbations 
in the CMBR, through the consideration of the Sachs- Wolf effect, describing the fractional 
variation of the temperature dT/T. Our procedure will not involve details of the radiation 
to matter era transition or of the matter era itself and so will only be applicable to the 
lower multipoles of the radiation; in particular we shall calculate the quadrupole moment 
to obtain the CMBR normalization. Here we briefly describe how we implement a rather 
standard method (see for example ref. |26[|) using the power function of if). (We have also 



calculated using another method, that of ref. |^TJ , which takes account of some details and 



have obtained similar results.) We track the development of the metric perturbation if), from 
the end of reheating, through the radiation era and into the matter era. 

The cosmic microwave background is described by a single number, the temperature 
T of the black body radiation. The small anisotropies in it from the direction dependent 
fluctuations, 5T(e), where e is a unit vector in the direction of observation, depend on the 
metric perturbation along the light path after emission at last scattering - the Sachs- Wolfe 
effect. This can be expressed as [p5l ,Effil 

^(e) = -^(x L5 ,r L5 ) (50) 
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where x^s = -J-e, -J^ being a good approximation to the distance of the last scattering sphere 
with H the present Hubble and tls the emission time. 1Z(x.ls, tls) is stochastic being 
given linearly by i/j(xls, tls) which is itself stochastic being expressed by the corresponding 
formula to Eq.(||). 

Using Eq.floTf) the angular correlation moments, Q, can be calculated in terms of IZk^TLs)- 
With ensemble averaging using Eq.(||) it can be shown, after some manipulation |p6 |, that 



C 



25 



dk 



V n {k) 



Ji 



' 2k ' 
a H 



1 2 



where V-jz(k) is the power function of the curvature perturbation: 

k 3 



V{k) 



2tt 



\n k \ 2 



(51) 



(52) 



To calculate low moments we only need TZkijis) for values of the order of aoH and 
at that time they are well outside the horizon, k <C aH{jLs) and so IZk ~ — (k to a good 
approximation. So we can replace 1Z by ( in Eqs.(^l|) and (|52| . Knowing (k at the end of 
inflation for an appropriate range of k enables the calculation of correlation moments using 
these equations. This can be done in principle by computing ( k for sufficient fc-values to 
cover the relevant range defined by the spherical Bessel functions. Alternatively we note 
from section |11 B| that ipk,Ck ^ kP" 1 ^ 2 at the beginning of reheating. The values of k/aH 
in the reheating period are extremely small for values of k appropriate to the CMBRF 
(and indeed even much larger values) and so the k 2 terms in the reheating equations are 
negligable. Thus the /c-dependence passes unaltered through the reheating and we use it to 
perform the ^-integration of Eq. ( |5T| ) to produce the results which now follow. 

We calculate the quadrupole moment defined, with Z = 2, as 



Q = T 



'21 + 1 



Air 



C, 



(53) 



Let = (k,out/ Ck,mi this ratio being independent of k for the relevant k and is found from the 
numerical solution of the reheat equations (Tables I and II). In terms of the /c-independent 
quantity Rc i C,k,ink 1 ^ 2 ~ p we find 



C-2 



i?dOMn|fc 1/2 " P (a H Q 



(54) 



p 25nJ z J2K ' 



From Eqs. [27|-[32] 

\Ck, n \k 1/2 - p (a H ) p+1 



2irp ( a H Q \ p+1 H, 



m Planck 



■\M(p)\ 



(55) 



(56) 



M(p) 



^{-p/2f 



sin((l/2-p)7r)r(p+l/2)' 



(57) 
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For g/m < 4500 and for the values of the decay constants given in Table I we find 



Q w !2{iK. 



(58) 



If ( had passed unchanged through reheating we would have found 



Q « UfiK 



(59) 



The observed value is Q m 18/jiK, and our theoretical result can readily be changed to 
agree with this by adjustment of the parameters y?o,m - the homogeneous inflaton field at 
the end of inflation - and to, the reheating potential parameter. Similar relatively small 
adjustments to these parameters (while scaling g,M, T x and T v proportionately to to) can 
give agreement with the observed Q for larger couplings such as those of Table II. 

Of course, while the evaluation above verifies that the parameters we have used are in 
the right region, the point is not the precise values of Q found but rather the amount of 
change induced in ( by the reheating which forces changes in the parameters of the theory. 



We have used a well-known simple theory that can give rise to parametric resonance of 
a scalar field, x> coupled to the scalar inflaton, </?, through a term g 2 (p 2 x 2 - We have inves- 
tigated preheating, followed by reheating mediated by decay terms T^ip 2 and T x y 2 in the 
Lagrangian which implement a smooth transition to a radiation era with a hydrodynamic 
fluid of relativistic particles. We have arranged a transition from an inflationary era ro 
the heating era that is also smooth. Besides non-linear effects of the extra field x we have 
included non-linear effects of the inflaton particle field ipi. The investigations centre on the 
three-curvature related parameter ^ for values of k relevant to the cosmic microwave back- 
ground fluctuations. While, for such k, (k is constant through the latter part of inflation and 
also through the radiation and subsequent eras as far as the generation of the fluctuations, 
the results presented show (k varying in the heating transition. Thus they emerge into into 
the radiation era with a value changed from that in inflation which implies a change in the 
fluctuations - or rather a change in the parameters of the theory from what they would have 
been if ( had been constant throughout. The change depends markedly on the magnitudes 
of g 2 , r v and T x . 

For given values of T v and T x there is a range of g 2 from small up to some value depending 
on the the magnitudes of T v and T x where parametric resonance occurs but where the 
change in ( is always by the same proportion. Though in the results presented this change 
is a diminution only of the order of 10% or so, it is interesting that it appears not to be 
due to the parametric resonance in x but only to non-linear effects of the inflaton particles 
produced - the energy density of the x particles never quite reaching that of (p particles. 
This change in ( should be no great surprise as when < ipi > is comparable with \<pl\, as it 
is during the latter stages of heating, theorems on the constancy of ( break down. However 
there is a range of g 2 , and T x where there is considerable parametric resonance but the 
view that ( should be unvarying throughout is a reasonable approximation. 

With larger g 2 then ( increases by a considerable proportion this being up to a factor of 
20 in the examples we have given. These cases are characterised by the x energy density, 



V. SUMMARY AND DISCUSSION 
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after the onset of parametric resonance, approaching and at times surpassing the <pi energy 
density. As can be seen from the figures , ( fluctuates markedly after 2 or 3 efolds into the 
heating transition. If we make the Ts larger so that the decay to radiation is quicker then 
the variability of the radiation era £ with g 2 increases. 

In the inflationary period the x field is restricted to within about the primordial value 
of quantum fluctuations. As seen from Figures 1 and 3 this persists for a few efolds into the 
heating transition before parametric resonance sets in and raises the energy density by about 
50 powers of ten in half an efold. Subsequently the x energy density never greatly exceeds 
that of the (pi, it seems that the direct interaction between these two fields is a strongly 
limiting factor. That is, the same interaction that gives rise to parametric resonance may 
also act to limit it because it contains inflaton particles as well as the classical inflaton field, 
(fo- This feature seems likely to apply to many other models with parametric resonance. 
So then the inflaton particles always play an important role in the evolution of ( and the 
revelation of their significance may be the most important result of the non-linear theory. 

The formalism we use has been explained at much greater length in a preceeding paper 
||10|| . In that paper also fairly strong ( variations were computed but the results given here 
supercede (see Section |ni|) those previous ones [fT0|JTTJ| . Also the present work extends the 
previous work by treating the completion of reheating into the radiation era where ( becomes 
a constant, different from that at the end of inflation, and thus the consequences for the 
cosmic microwave background fluctuations become evident. Our model extends smoothly 
from the primordial quantum excitations to the creation of the fluctuations. 

The main feature of our formalism is that as the number density of the original quantum 
fields increases a smooth transition to where they can be treated as classical stochastic occurs 
T8fl ; the only manageable treatment of the non-linear terms in the equations of motion is to 



take ensemble averages. This gives rise to loop imtegral factors, such as that of Eq.(^), in 
the equations of motion. The evaluation of these involves both an infra-red, low k, cut-off 
and an ultra-violet, high k, cut-off. This has been discussed in detail in Appendix B of 



ref. 10] where the low k cut-off was regarded as a contribution to the renormalization of 



the homogeneous part of the field. The high k cut-off is fixed to correspond to a wavelenth 
H -1 , the inverse Hubble length at the beginning of heating. Such a cut-off has been used 



in the work of other authors for example in ref. [ 12j . Though the only obvious approximate 
magnitude it is a major source of uncertainty in our work, as any cut-off very much less 
or very much greater would significantly affect our results. A very much smaller cut-off of 
course corrsponds cloosely to the neglect of all non-linear effects. 

Where, on the contrary, we have used a linear approximation is in working in the equa- 
tions of motion to first order in the metric perturbation if). The measure of validity we 
have taken is, in accord with the metric Eq.(|l|), to compare 2\/< if> 2 > with unity fllPf . In 
cases where ( decreases this corresponds throughout the heating process to a less than 1% 
perturbation to the metric and the first order perturbation theory is well justified. But in 
cases where ( increases it can reach up to nearly a 10% perturbation in the latter stages of 
reheating. While we do not think that this gives doubt that these are cases where ( increases 
markedly nevertheless the quantitative aspect would benefit from a more refined treatment. 

We find it difficult to see any direct observational test of the changes in the heating 
transition. In the body of the paper we have used that the shape of the power spectrum, in 
our theory a power law, emerges unchanged from the transition and so only depends on the 
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inflationary theory. The parameters of nearly all envisaged theories can be changed to follow 
the magnification (or diminution) of ( and adjust to the COBE or other normalization. The 
dependence on the details of preheating and reheating makes this adjustment uncertain in 
the present state of knowledge. The best chance may lie in a comparative calculation of 



gravitational wave production in the hope that this can be observed [28,29 
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TABLES 



TABLE I. The change in ( in the inflation to radiation era transition, shown for smaller values 
of the X) V coupling, g. Also shown are the increase, n(efolds), in ln(a) in the transition and the 
ratio of densities of x an d V?i at the maximum density of x- The decay constants are given by 
T x /m = T v /m = 5 x 10" 5 . 



g/m 


n(efolds) 


log (denx / dentpi ) 


Cout/ Cin 


1000 


7.9 


-39 


0.88 


3000 


7.6 


-3.9 


0.87 


4000 


7.9 


-14 


0.88 



TABLE II. The increase of C m the inflation to radiation era transition, shown for some larger 
values of the x, if coupling, g. Also shown are the increase, n(efolds), in ln(a) in the transition 
and the ratio of densities of x and ipi at the maximum density of x- The decay constants are given 
by T x /m = T^/m = 5 x 10~ 5 . 



g/m 


n(efolds) 


log(denx/denipi ) 


Cout/ Cin 


4700 


7.9 


-.9 


5.6 


4900 


7.9 


-.9 


10.2 


5000 


7.6 


-.5 


19.1 


7000 


7.5 


+.7 


6.3 


12000 


7.9 


-.9 


10.3 
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FIGURES 




FIG. 1. For g/m = 3000, logarithmic plot of energy densities of x (A)and ip\ (B) versus the 
number of e- folds of expansion in the reheat era. In all figures densities are in units m 2 x m 2 Planck , 
where m = KT 7 m„ ancfc . 




3 4 5 

n(cfolds) 



FIG. 2. For g/m 
reheat era. 



3000, plot of lg(C/Cm) versus the number of e- folds of expansion in the 
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n(efolds) 

FIG. 3. For g/m = 3000, logarithmic plot of energy densities of ipi (A), ipo (B) and po versus 
the number of e-folds of expansion in the reheat era. 




4 5 
n(efolds) 

FIG. 4. For g/m = 4700, logarithmic plot of energy densities of x (A)and 991 (B) versus the 
number of e-folds of expansion in the reheat era. 
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n(efolds) 

FIG. 5. For g/m = 4700, plot of lg(C/Cm) versus the number of e- folds of expansion in the 
reheat era. 




012345678 

n(cfolds) 

FIG. 6. For g/m = 4700, logarithmic plot of energy densities of ip\ (A), ipo (B) and po versus 
the number of e-folds of expansion in the reheat era. 
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